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We define an activity dependent branching ratio that allows comparison of different time series 
Xt- The branching ratio b x is defined as b x = E[£ x /x]. The random variable £ x is the value of 
the next signal given that the previous one is equal to x, so £ x = {Xt+i\X t = x}. If b x > 1, the 
process is on average supercritical when the signal is equal to x, while if b x < 1, it is subcritical. 
For stock prices we find b x = 1 within statistical uncertainty, for all x, consistent with an "efficient 
market hypothesis". For stock volumes, solar X-ray flux intensities, and the Bak-Tang-Wiesenfeld 
(BTW) sandpile model, b x is supercritical for small values of activity and subcritical for the largest 
ones, indicating a tendency to return to a typical value. For stock volumes this tendency has an 
approximate power law behavior. For solar X-ray flux and the BTW model, there is a broad regime 
of activity where b x — 1, which we interpret as an indicator of critical behavior. This is true 
despite different underlying probability distributions for Xt, and for £ x . For the BTW model the 
distribution of £ x is Gaussian, for x sufficiently larger than one, and its variance grows linearly with 
x. Hence, the activity in the BTW model obeys a central limit theorem when sampling over past 
histories. The broad region of activity where b x is close to one disappears once bulk dissipation is 
introduced in the BTW model - supporting our hypothesis that it is an indicator of criticality. 

PACS numbers: 89.65.Gh, 89.75.Da, 89.75.Fb, 05.45.Tp, 05.65.+b 



I. INTRODUCTION 

Detailed forecasting in complex systems is often dif- 
ficult if not impossible. Nonlinear processes as well as 
long range spatial and/or temporal correlations can ren- 
der a direct, reductionist approach futile. Furthermore, 
in many cases of interest, controlled laboratory experi- 
ments are unfeasible. For instance, stock market data or 
solar X-ray flux can only be obtained under specific con- 
ditions set by the system itself, and observing the time 
series under various controlled conditions is not possible. 

Testing the efficient market hypothesis (EMH) 
presents a clear example of this difficulty. Roughly speak- 
ing, the EMH states that asset prices are inherently un- 
predictable [H 0,O| > or that the market is hard to beat 
There are many flavors of the EMH. The weak EMH 
states that the market is efficient if agents only have in- 
formation about the time series of market prices. The 
strong EMH, on the other hand, states that the market 
is efficient when agents have access to all relevant infor- 
mation that could affect prices; this includes e.g. insider 
trading. 

Mathematically the weak EMH can be formulated in 
terms of a martingale property for the time series of 
prices p|. In its simplest form, a stochastic variable, X t , 
is said to be a martingale if the expectation of its next 
value given its entire past is equal to its current value, 
or E[X t +i\X t ...Xo] = X t for all t. Empirically, it is not 
possible to obtain this expectation value directly from 
real world time series. In fact, the existence of EMH in 
any of its forms is highly disputed [1, 0] ■ 

For a Markov process, the value that X t +\ takes only 
depends on the previous value X t . Indeed a necessary but 



not sufficient requirement for any stochastic process to be 
a Martingale is that i?[A t+ i|X t ] = X t . One example is 
the critical Galton- Watson (GW) branching process 
Starting with a single node, Xq = 1, each node inde- 
pendently produces a number of offspring that is Poisson 
distributed with mean b. Here b is called the branch- 
ing ratio. If b = 1 the process is critical and is also a 
Martingale, since E[X t+1 \X t ...X ] = E[X t+1 \X t ] = bX t . 

If the underlying probability distribution used to eval- 
uate the expectation value E[- ■ •] is not known, or if the 
nodes interact with each other in generating offspring, 
one can still empirically measure an activity dependent 
branching ratio as b x — E[^ x /x\. Here the random vari- 
able £ x is the value of the next signal given that the previ- 
ous one is equal to x, or £ x = {X t +i\X t = x}. We inter- 
pret t; x as the set of outcomes of an interacting branching 
process with a current population x. Empirically the ex- 
pectation value E[£ x /x] is an average over all times t 
when X t = x. If b x = 1 the process is on average critical 
when the activity is equal to x. If b x > 1, it is supercrit- 
ical and if b x < 1 it is subcritical. Note that we are not 
making any assumptions that the processes we measure 
are, in fact, Markovian. The measured branching ratio 
b x is an average over all observed histories leading to a 
population of size x. 

We use the activity dependent branching ratio b x to 
compare and contrast time series from stock markets, a 
physical system (solar X-ray flux) and a model (the Bak- 
Tang-Wiesenfeld (BTW) sandpile). Previously, time se- 
ries of activity in the BTW sandpile have been compared 
in detail with that of solar flux - finding a number of 
similarities 0. Our analysis finds similarities as well as 
significant differences in these two systems. 
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We find that b x is statistically indistinguishable from 
unity for time series of stock prices as well as for the Dow 
Jones industrial average, consistent with the weak EMH. 
On the other hand, stock volumes, X-ray flux and activ- 
ity in the BTW model all show roughly similar behavior: 
b x decreases from a supercritical value at low levels of ac- 
tivity to a subcritical one at large values. This indicates 
a general tendency for the activity to return to a charac- 
teristic value, which is not present for stock prices. For 
stock volumes, the branching ratio has a relatively strong 
dependence on activity, by ~ (V/{V))~ a with a ~ 0.69. 
Solar X-ray flux and activity in the BTW model both 
show a broad range of activity where the branching ratio 
is close to one. This broad range increases with the sys- 
tem size for the BTW model and disappears once bulk 
dissipation is introduced, suggesting that it is an indica- 
tor of criticality. 

We also compare and contrast the probability distribu- 
tions P(^ x /x), at particular values of activity x in these 
systems. For the BTW model, this distribution is well 
described by a Gaussian for x sufficiently larger than 1. 
On the other hand, for both stock volumes and flux inten- 
sities, P(£ x /x) is broader than Gaussian. The marginal 
distribution of flux intensities P(I) is well described by a 
power law, while for the BTW model the distribution of 
activity P{n) is approximately exponential with a corre- 
lation length that grows with system size. 

In Section II, we present results for the branching ratios 
determined from analyses of time series for stock prices 
and for stock volumes considering four different stocks as 
well as the daily Dow Jones average. Section III presents 
results for solar X-ray flux data. Section IV contrasts and 
compares results for the canonical BTW model with vari- 
ants (a) including bulk dissipation - making the model 
subcritical, and (b) having periodic boundary conditions, 
which leads to non-ergodic behavior. We also discuss how 
our activity dependent branching ratio differs from the 
average branching ratio for avalanches in self-organized 
critical (SOC) systems previously discussed in the lit- 
erature (See e.g. Refs. (HI)- Section V contains a 
discussion and summary of the main results. 



symbol (■ ■ •) indicates an average over the observation 
time. For all price time series studied b v = 1 within 
statistical uncertainty. These data are binned such that 
there are at least 500 points in each bin, and error bars 
indicate one standard deviation. 
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FIG. 1: (color on-line) The branching ratio b p vs. price for 
different stocks. The x-axis has been scaled by the mean 
price. All data shown have a resolution of one minute, where 
the price used is that at the start of the minute. Error bars in 
this and subsequent figures indicate one standard deviation. 
The data for qqqq was taken over the period 09:30 23/05/08 

- 13:37 20/06/08, the data for csco is from 09:30 23/05/08 

- 13:37 20/06/08, the data for aapl is from 9:30 27/05/08 - 
14:08 23/06/08, the data for F is from 09:30 27/05/08 - 14:03 
23/06/08. These results are consistent with the weak EMH. 
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II. STOCK MARKET 

Does knowledge of the current price or volume of trade 
for a particular stock or market index allow one to make 
predictions about the next value? For prices the answer 
is no, while for volume of trade the answer is yes, as 
discussed next. 

We analyze one minute resolution data for four differ- 
ent stocks from [l2j for intervals of 28 days. We also 
examine one day resolution data for the Dow Jones over 
80 years from [l|. Both price £ P = {P(t + l)\P(t) = P}, 
and volume £y = {V(t + l)\V(t) = V} are studied. 

Fig. Q] shows the activity dependent branching ratio 
bp = E[£ p /P] vs. P/{P) for four different stocks. The 
same quantity for the Dow Jones is shown in Fig. [21 The 



FIG. 2: (color on-line) Same as Fig. 1 for the Dow Jones in- 
dustrial average. The data is from 01/10/1928 - 23/05/2008, 
and has a resolution of 1 day using the opening price. The 
behavior is also consistent with the weak EMH. 

The activity dependent branching ratio for volume, by, 
has a strong dependence on volume as shown in Fig. [3] 
For small values the stocks behave like a supercritical 
branching process, while for large values they are sub- 
critical. Hence the volume has a tendency to return to 
roughly its mean value. In fact, by has an approximate 
power law dependence on V, with by ~ (V/ (V))~ a . The 
exponent a ~ 0.69 for three of the stocks shown in Fig. [3] 
but appears to be smaller (or nonexistent) for the Ap- 
ple stock (aapl), which also has a more limited variation 
in volume, precluding any firm conclusion about scaling. 
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We have also analyzed this quantity for different time 
windows (data not shown) and the results for the individ- 
ual stocks do not vary in any substantial way. The data 
for the Dow Jones volumes (not shown) are too noisy to 
draw definite conclusions. 
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FIG. 3: (color on-line) Activity dependent branching ratios 
for stock volumes during the same time period as in Fig. [I] A 
line with slope m = —0.69 is included as a guide for the eye. 
The behavior of the Apple stock (aapl) appears to deviate 
from the other three. 

For a given level of activity, the probability distribution 
P(£v/V), also differs significantly from that for prices 
P{£p/P). For prices we find results similar to that found 
in Ref. [14], who analyzed price changes for different time 
increments over all prices. In our case the data (not 
shown) are much more noisy since we fix both the initial 
value of price, as well as the time interval (one minute). 
For stock volumes and solar intensities, the cumulative 
distributions P(£v/V > b) at a given V and P(£i/I > b) 
at a given /, are both broader than Gaussian as shown 
in Fig. II 



III. SOLAR X-RAY FLUX 

Solar flares are bursts of radiation that occur in the 
solar corona. These bursts can reach sufficiently high 
energies to pose a risk to astronauts, spacecraft, or 
airplanes following polar routes. In addition they ex- 
hibit a number of emp irical features associated with 
SOC @, GJ, GJ, HI EE El- For instance, the distri- 
bution of event durations and quiet times is a power-law 
for both solar flux and the BTW sandpile [9j, once physi- 
cally relevant detection thresholds are taken into account 
to compare these time series on an equal basis. Our anal- 
ysis shows that the dependence of the branching ratio on 
activity are similar in the two cases although the under- 
lying probability distributions for activity (X t ) and for 
subsequent conditioned activities £ x are different. 

We examine time series in the 1 — 8A range ob- 
tained from the Geostationary Operational Environmen- 
tal Satellites (GOES) Satellites at the "Space Physics 
Interactive Data Resource" [20]. A time series of five 
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FIG. 4: (color on-line) Comparison of cumulative distribution 
function (CDF) P(£x/x > b) for stock volumes and solar in- 
tensities at a given value of activity x. For stock volumes, 
Ford (F) was used, and the CDFs calculated at x = V/ (V) = 
0.2 ± 0.05 (supercritical region) and x = V/ (V) = 2 ± 1 (sub- 
critical). For solar intensities the entire time series, "AH", was 
used, and the CDFs calculated at x = 1/(1) = 0.1 ±0.005 (su- 
percritical) and x = 1 1 '(/} — 100 ± 10 (subcritical) . The Ford 
CDF at x — 0.2 ± 0.05 is compared with a Gaussian having 
the same mean and variance. The distributions are broader 
than Gaussian in all cases. 
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FIG. 5: (color on-line) Activity dependent branching ratios 
for solar flux intensity at solar minimum, solar maximum, 
and for the entire data set, "All", as defined in the text. This 
shows a weak tendency to return to a typical value although 
bi ~ 1 for a broad range of intensities, /. This behavior is 
comparable to that shown for the SOC BTW model in Fig. [7] 



minute intervals from 01/01/1986 to 30/04/2008 was ob- 
tained from GOES satellites 5-12. When data from mul- 
tiple satellites was available the average of the available 
data was used. The time series spans approximately two 
solar cycles. Data from 10 6 — 3 x 10 6 minutes were taken 
to correspond to a solar maximum, and the portion from 
4 x 10 6 — 6 x 10 6 minutes to a solar minimum. 

Only values of / > Iq > 3 x 10~ 8 W/m 2 were used when 
computing statistics. This value is close to the detection 
threshold of the satellites. In order to make a compari- 
son with the BTW data the solar X-ray intensities were 
divided by Iq. This transforms the minimum possible 
value in both data sets to one. 

The behaviour of 6/ is qualitatively similar to by, as 
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FIG. 6: (color on-line) The probability distribution function 
of solar x-ray intensities, P(I). The straight line indicates a 
power law fit with exponent 2.3 for the "All" time series. This 
is different from the approximately exponential behavior seen 
in the SOC BTW model as shown in Fig. U 



the system pi I22] . 

A. Self-organized Critical BTW model 
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shown in Fig. [5J The branching ratio decays from su- 
percritical to subcritical as I increases. However, unlike 
by, bi is close to one over a broad range of intensities, so 
the tendency to return to a characteristic value is much 
weaker for flux intensities than for stock volumes. This 
broad range where b « 1 is also a property of the SOC 
BTW model as shown in the next section. It disappears 
once bulk dissipation is introduced into the BTW sand- 
pile; hence we interpret this broad range as a signal of 
critical behavior. 

Despite this close similarity, the BTW model and solar 
activity drastically differ with respect to the distribution 
of activity, P(X t ), and the distribution of (£ x /%) a t a 
given activity. As shown in Fig. [6j the probability dis- 
tribution function for flux intensities P(I) is broad with 
a tail that is consistent with a power law with exponent 
~ 2.3. On the other hand, the probability distribution 
of activity P(n) in the BTW sandpile, shown in Fig. 
is close to, but not exactly, exponential. In addition the 
cumulative distribution function P(£//I > b) at a given 
level of activity / is broad as indicated in Fig. 2] This 
contrasts with the BTW model, where the distribution 
P(£ n /n) at a given level of activity, n, is Gaussian (see 
Fig. IS). 

IV. BTW SANDPILE MODELS 

The BTW sandpile model [2l[ is the paradigmatic ex- 
ample of SOC. SOC describes slowly driven, dissipative 
systems that reach a critical state without fine tuning 
parameters. The BTW model depicts a system that is 
externally driven to a local dissipative instability where- 
upon it "topples" . This toppling can induce further 
topplings, which can lead to cascades of activity prop- 
agating through the system. These cascades are called 
avalanches. In the steady state, the BTW model reaches 
a stationary state where the distribution of avalanche 
sizes is broad with no natural scale other than the size of 



FIG. 7: (color on-line) Activity dependent branching ratio for 
the SOC BTW model for different system sizes, L. Error bars 
are smaller than symbol size. As the system size increases, 
the region where b„ m 1 broadens. The entries beginning 
with D denote the dissipative BTW model. The dissipative 
"D-BTW" model does not show a broad region where b„ ~ 1. 

The SOC BTW sandpile model is composed of an L x L 
lattice with open boundary conditions, where each site 
is assigned a height z. The height of a stable site is an 
integer between zero and three. A site with a value z > 3 
becomes unstable and topples by adding a grain to each 
of its four nearest neighbors, thus decreasing its height by 
four. If a boundary site topples, it throws some grain(s) 
out of the system. Initially the sandpile is empty and 
z = for all lattice sites. The system is driven by adding 
a grain to a randomly chosen site. Then all unstable sites 
are updated in parallel, and the time unit is increased 
by one. This continues until all sites are stable. Then 
another grain is added and the process is repeated ad 
infinitum. We start collecting statistics after the average 
number of grains in the pile becomes stationary. A time 
step corresponds to one parallel update of all lattice sites, 
or to the addition of a single grain, whichever is the case. 

At every time step, i, we record the number of toppling 
sites, n t . We define an activity dependent branching ra- 
tio b n = E[nt+i/n\n t = n], as the fraction of sites that 
topple in a time step immediately following one where n 
sites topple. We have numerically simulated the BTW 
model on lattice sizes ranging from L = 100 to L — 6000 
to study finite size effects. 

Fig. [7| shows b n vs. n/{n) for various system sizes for 
the SOC BTW model. The behavior is qualitatively sim- 
ilar to that for solar flare intensities, and for stock market 
volumes. The range of n where b n is close to one in- 
creases with system size L. The subcritical region occurs 
for large n as a result of dissipation at the boundaries, 
which limits the maximum size of n. This means that 
after a large toppling event the system is more likely to 
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undergo a smaller one and b n < 1. We attempted a fi- 
nite size scaling analysis, which did not give compelling 
results. This is consistent with previous results indicat- 
ing that the SOC BTW model does not obey finite size 
scaling [HHl. 

Fig. [8] shows the probability distribution P(n), which 
has an approximately exponential decay. This differs 
from the comparable result shown in Fig. [5] for the solar 
X-ray intensity. Hence, the behaviour exhibited by b x is 
robust for systems that have markedly different distri- 
butions for activity. As shown in Fig. [5J we attempted 
a finite size data collapse of the distribution of activity, 
but this collapse shows systematic deviations. However, 
it is clear that for the distributions the correlation length 
increases with system size, leading to a broadening distri- 
bution of activity in the large L limit for the SOC BTW 
model. 
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FIG. 8: (color on-line) The probability distribution for activ- 
ity in the SOC BTW model, P(n). Error bars are smaller 
than symbol size. The decay is approximately exponential. 
The best data collapse of the tail of the distribution is ob- 
tained by rescaling with L 1 ' 2 . Since the SOC BTW model 
does not exhibit finite size scaling this rescaling is only for 
the purpose of plotting all the data together. 

We examined the probability distribution, P(£„/n) for 
various values of n. For n sufficiently larger than one, the 
distribution is well-described by a Gaussian. Fig.[9]shows 
the distribution function for three values of n in a system 
of size L = 500, one in the supercritical regime and two 
in the subcritical one. Fig. ITD1 shows that the variance of 
the random variable increases linearly with n. Gaus- 
sian behavior with a variance that grows linearly with n 
indicates that the activity in the BTW model obeys a 
central limit theorem: when sampling over prior histo- 
ries, for each n the activity is the sum of n independent 
random processes. 

B. Dissipative BTW Model 

We analyze a BTW model that includes bulk dissipa- 
tion to test how criticality affects the activity dependent 
branching ratio. The model is similar to the SOC BTW 
model except it also includes bulk dissipation (2f| . When 
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FIG. 9: (color on-line) Probability distribution function 
P(£ n /n) for the SOC BTW model. Error bars are smaller 
than symbol size. The distribution is shown for n = 
15, 150, 450, with L = 500. The system for n = 15 lies in 
the supercritical region, while for n = 150 and n = 450 it 
lies in the subcritical region. In all cases the distributions are 
indistinguishable from Gaussian. 
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FIG. 10: (color on-line) The variance of the random variable 
£n, cr 2 (C")/ n vs. n for three different system sizes. The ratio 
goes to a constant for large n. Hence the variance of f„ grows 
as n. 

a site topples all its grains are removed from the system 
with probability pd, and with probability 1 — pd the nor- 
mal toppling rule applies. Fig. [7] compares the branching 
ratio for the dissipative BTW model with pd = 10 -2 to 
the SOC version. It shows that the broad region where b n 
is approximately equal to one disappears once dissipation 
is introduced. 



C. BTW Model with Periodic Boundary 
Conditions 

We also studied the BTW model with periodic bound- 
ary conditions, so that no grains are ever thrown out of 
the system. As grains are added, an infinite avalanche 
eventually occurs. We only examine statistics of the 
infinite avalanche. This corresponds to a fixed energy 
sandpile, which have been previously studied in [26l . [271 . 
|28| . [29(| . In our analysis, an avalanche that lasts more 
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than 9 x 10 7 parallel update steps is considered to be 
infinite, and we only collect statistics during the infinite 
avalanche. 

Fig. [TT] shows the time series n t for four infinite 
avalanches during 6000 time steps. The figure shows 
that nt is periodic and not ergodic, as was previously 
noted in e.g. Ref. 29]. For each realization of the infi- 
nite avalanche the dynamic range of n is small compared 
to SOC BTW model. Moreover, the system is sensitive 
to initial conditions. We tested this sensitivity to initial 
conditions by starting the lattice empty, or by randomly 
initializing each site to a value of or 1. The initial con- 
ditions affects both the period, and the amplitude of os- 
cillations of the infinite sized avalanche. Similar results 
were also obtained by simply keeping the same initial 
conditions and changing the seed of the random number 
generator used. Changing the seed alone was enough to 
similarly affect the period and amplitude of oscillations. 
All these results imply that the BTW model with peri- 
odic boundary conditions is not ergodic and cannot be 
compared to the other systems studied here in a mean- 
ingful way. 




t(parallel updates) 

FIG. 11: (color on-line) Time series of activity n t for different 
initial conditions for the BTW model with periodic boundary 
conditions. This displays different oscillatory behaviour for 
different initial conditions. In each case a system of size L — 
500 was used. 



D. Previous definitions of the branching ratio 

The activity dependent branching ratio defined in this 
paper differs from the average branching ratio measured 
in [Tol . [Tl| . The previously defined ratio was not con- 
ditioned on activity but rather defined as the average 
activity, over all avalanches, resulting from a single top- 
pling. Indeed it was shown that this average branching 
ratio b = 1 — l/(s), where (s) is the average avalanche 
size. Hence b is not an independent quantity, and is al- 
ways, by definition, less than or equal to one, as long as 



the average avalanche size is finite. Our activity depen- 
dent branching ratio is not restricted to situations where 
an avalanche can be well-defined or one can identify in- 
dividual sites for activity. In addition it gives an overall 
picture for how the system behaves at different levels of 
activity, unlike the average in Ref. [l(| EH , which sums 
over all observed levels of activity. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we introduced an activity dependent 
branching ratio, b x , and use it to analyze different time 
series, X t , arising in physical, economic, and model sys- 
tems. 

We found that stock prices have a branching ratio in- 
distinguishable from unity over for all observed prices. 
This observation is consistent with the weak efficient 
market hypothesis. Conversely, stock volume, solar X- 
ray flux, and the self-organized critical BTW model ex- 
hibit supercritical branching ratios for small levels of ac- 
tivity and subcritical ratios for large ones. This indicates 
a tendency for these systems to return to a characteris- 
tic value. This tendency is most pronounced for stock 
volumes which show a trend consistent with power law 
with exponent ~ 0.69, for three out of four of the stocks 
examined. It is not yet clear what separates the Apple 
stock in our analysis from the other three, or what the 
source of the apparent scaling is. Solar X-ray flux, and 
the BTW model both show a broad region where the 
activity dependent branching ratio b x « 1. When bulk 
dissipation is introduced into the BTW model this broad 
region disappears, supporting our hypothesis that this is 
a signature of criticality. 

The BTW model and solar X-ray flux show this simi- 
larity despite having different underlying probability dis- 
tributions for X t . For solar X-ray flux the distribution 
of flux intensities is consistent with a power law with ex- 
ponent ~ 2.3, while for the self-organized critical BTW 
model the distribution of activity P{n) has an approxi- 
mately exponential decay, with a correlation length that 
grows with system size. 

We also found that the variance in activity <7 2 (£„) 
scales linearly with n for the BTW model, and the dis- 
tribution of subsequent activity is Gaussian at a fixed n. 
This indicates that the BTW model obeys a central limit 
theorem when sampling over past histories. It remains 
to be seen if this last result can be derived theoretically. 
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